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Abstract 

We consider the integral light polarization from optically 
thick accretion disks. Basic mechanism is the multiple light 
scattering on free electrons (Milne's problem) in magnetized 
atmosphere. The Faraday rotation of the polarization plane 
changes both the value of integral polarization degree p and 
the position angle \- Besides, the characteristic spectra of 
these values appear. We are testing the known relation be- 
tween magnetic field of black hole at the horizon Bbh and its 
mass Mbh-i and the usual power-law distribution inside the 
accretion disk. The formulae for p(X) and xW depend on a 
number of parameters describing the particular dependence 
of magnetic field in accretion disk (the index of power-law 
distribution, the spin of the black hole, etc.). Comparison of 
our theoretical values of p and x with observed polarization 
can help us to choice more realistic values of parameters if the 
accretion disk mechanism gives the main contribution to the 
observed integral polarization. The main content is connected 
with estimation of validity of the relation between Bbh and 
Mbh- We found for the AGN NGC 4258 that such procedure 
does not confirm the mentioned correlation between magnetic 
field and mass of black hole. 

Keywords: polarization - magnetic fields - accretion disks 
- supermassive black holes. 

1 Introduction 

Recently Zhang et al. (2005) used the observed ^L(5100/l) — 
Mbh correlation to probe magnetic field of black holes har- 
bored in AGNs. Their model uses the assumption that the 
accretion disk is supplied by the energy injection due to 
the magnetic coupling (MC) process, and the gravitational 
dissipation due to accretion. They have modeled the rela- 
tion vL(51Q0A) — Mbh as a function of the spin parameter 
a I Mbh and the magnetic field Bbh- Using the observations 
of i/L (5100 A) for 143 AGNs, they obtain the correlation be- 
tween Mbh and Bbh of the following form: 



log B BH (G) = 9.26 - 0.81 log 
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According to this formula, for Mbh/Mq = 10 , 10 , and 
10 8 the magnetic field Bbh is equal to 1 G, 100 G, and 630 
G, respectively. Note that Zhang et al. (2005) have used 
the theory of magnetic coupling (see Wang et al. 2002, 2003 
and Page & Thorne 1974), which supposed that magnetic 
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field is poloidal with the dominant vertical component to the 
accretion disk, i.e. in Eq.(l) is assumed that Bbh = Bu. 

Relation (1) differs significantly from the relationship based 
on the standard accretion disk theory of the inner region (see 
Shakura and Sunyaev 1973; Novikov and Thorne 1973): 



B(G) = 10 8 
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where R g = 2GMbh /c 2 is gravitational radius of a black hole. 
Eq. (2) is based on the equipartition between the radiation 
pressure and magnetic pressure at the high accretion rate. 

For Shakura-Sunyaev zone (b), where the equipartition be- 
tween the gas pressure and magnetic pressure occurs, the 
magnetic field strength is determined by another relation: 
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where a is the Shakura-Sunyaev viscosity parameter, M and 
M^dd are the ordinary and Eddington accretion rates, respec- 
tively. 

The typical value of the magnetic field in the optical radi- 
ation region rj 10 3 i? g and for a = 0.1 is 
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Narayan and Yi (1995) proposed advection-dominated ac- 
cretion flow (ADAF) model. Their model presents the hot 
advection-disk solution that does include advection and then 
dominates the ion energy transfer. In this case a global mag- 
netic field can be presented as 
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The basic goal of our paper is to calculate the spectrum of 
linear polarization that can be used to determine the magnetic 
field strength in the region where the optical linear polariza- 
tion originates. We calculate the degree and position angle 
of the polarization of radiation scattered in accretion flows 
(standard accretion disk, ADAF model, etc.). We take into 
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account the Faraday rotation of the polarization plane (see 
Dolginov et al. 1995; Gnedin and Silant'ev 1997; Gnedin et 
al. 2006). The comparison with the observational data can 
allow us to test both the various models of accretion onto a su- 
permassive black hole and also the correlation between black 
hole mass and magnetic fields in AGNs and QSOs. Below we 
consider the case of relation (1) in more detail. 

2 Polarization from magnetized ac- 
cretion disks 

An accretion disk is one of the basic elements of the structure 
of accretion flows near supermassive black holes. Since there 
is no axial symmetry with respect to the line of sight, the 
integral, scattered by free electrons radiation, emerging from 
an accretion disk is polarized. The angle of Faraday rotation 
of a light beam passing through the distance I is described 
by the expression (see, for example, Dolginov et al. 1995; 
Gnedin and Silant'ev 1997): 



* = -StcosO ~ 0.4 ( 
2 \ 
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where r is Tomson's optical depth of I, A is the wavelength 
of the radiation and 9 is the angle between the magnetic field 
B and line of sight n. 
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is the Faraday dimensionless depolarization parameter. Here 
w = 2ttv is angular frequency, and lob = eB/m e c is cyclotron 
frequency of an electron in a magnetic field, r e = e 2 /m e c 2 is 
classic electron radius. 

The Faraday rotation of the polarization plane in a mag- 
netized accretion disk changes both the integral degree of po- 
larization p and the position angle \. The spectra of these 
values acquire the characteristic forms. 

Considering optically thick accretion disks we use the solu- 
tion of the Milne problem. Usually one considers the Milne 
problem corresponding to sources of thermal radiation being 
located far from optically thick plane atmosphere. The exact 
numerical solution of this problem for magnetic field Bp di- 
rected along the normal N to the atmosphere are presented 
in Agol and Blacs (1996) and in Shternin et al. (2003). 

For estimations of magnetic fields in accretion disks we use 
simple asymptotical formulae (see Silant'ev 2002 and 2005; 
Gnedin et al. 2006, and Silant'ev et al. 2009). These formulae 
are valid for arbitrary directions of magnetic field. Firstly we 
give the asymptotic formulae for Stokes parameters /, Q and 
U for general case of magnetized and absorbing plasma: 
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Here, F is the radiation flux emerging from the surface of 
accretion disk, fi = cosz, i is the angle of inclination of ac- 
cretion disk, i.e the angle between the outer normal N to 
the disk and the direction n to the observer. The value 
J(/i) = = 0) describes the angular distribution of 

escaping radiation intensity, J\ is the first moment of J(/i), 
i.e. the integral from /xJ(/x) over interval (0, 1). Parameter q is 
the degree of absorption. Dimensionless parameter C arises in 
turbulent magnetized plasma (see Silant'ev 2005), and char- 
acterizes the new effect - additional extinction of parameters 
Q and U due to incoherent Faraday rotations in turbulent ed- 
dies. Below we present the apparent formula for this param- 
eter. The value k is the root of corresponding characteristic 
equation. For conservative atmosphere (q = 0) this parameter 
is equal to zero. The parameter g characterizes the maximum 
degree of polarization p ma x = (1 — #)/(l + g), if we take into 
account only last scattering of non-polarized radiation be- 
fore escape the surface, neglecting the existence of magnetic 
field. So. for conservative atmosphere g — 0.83255. This 
corresponds to p max = 9.14%. Remember that the Sobolev - 
Chandrasckhar value is equal to 11.71%. 

The asymptotic formulae (8) suppose that the depolariza- 
tion parameter 5 is large. In this case the terms with parame- 
ters Q and U in full system of transfer equations for /, Q and 
U become very small <~ 1/(5, and they are negligible in the 
equation for intensity /. As a result, the radiation intensity 
obeys the separate transfer equation with the Rayleigh phase 
function (see, in more detail, Silant'ev 1994). The first equa- 
tion in system (8) presents the exact solution of this equation 
for the Milne problem. Formulae (8) were obtained in this 
way. 

The numerical values of k, g, J\ and J(p) for various values 
q are presented in Silant'ev (2002). So, for example, for q = 
0.1 we have k = 0.5232, g = 0.77129, p max = 12.91%, and for 
q = 0.2 the corresponding values are k = 0.70483, g = 0.70405 
and Pmax — 17.37%. It is known, that the polarization de- 
gree in Milne's problem in absorbing atmosphere is greater 
than in conservative one. So, for atmospheres with q = 0.1 
and q = 0.2 the p mtM; -values are equal to 20.4% and 28.7%, 
respectively (see Silant'ev 1980). In these cases the relative 
contribution of parameters Q and U into the polarization it- 
self is greater than for the case q — (where this contribution 
is ~ 20%). The existence of true light absorption can ex- 
plain the observed polarization degrees greater than limiting 
Sobolev - Chandrasckhar value 11.71% for conservative at- 
mosphere. Clearly, the parameter q depends on wavelength. 
This gives an additional mechanism to explain the wavelength 
dependence of observed polarization degree p(X) and position 
angle x(A). 

Usually one observes the axially symmetric accretion disks 
as a whole. The observed integral Stokes parameters 
(Q(n, B)) and (U(n, B)) are described by the azimuthal av- 
eraged formulae 
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Here Q(/x,0) - - g)(l - A* 2 )/2ttJi(1 + g)(l - M 

is the Stokes parameter of scattered non-polarized radia- 
tion (U(fi, 0)) = 0) in the classic Milne problem (see Chan- 
drasekhar 1950). The minus sign denotes that the electric 
field oscillations are parallel to accretion disk plane (it means 
that the position angle x(/x, 0) = 0). Further we denote 
p(/i,0) = |QOi,0)|. 

For estimations in conservative atmospheres [q = 0) in- 
stead of p(fx, 0) one can use pr{p) - the known value of 
polarization for non-magnetized atmosphere with Thomson 
(Rayleigh) scattering, which is presented in Chandrasekhar's 
book (1950). In this case we can also use small values of 
parameter 5, because the limiting value for 5 -C 1 is pr (/•*)• 

Dimensionless parameters a and b are connected with the 
parallel (along the normal N to the disk) B|| and perpendic- 
ular components of magnetic field B = By + B^: 

a = (1 - q)6\i cosi ~ 0.8(1 - q)\ 2 {nm)B\\ (G) p, 

b = (1 - q)8j_ sini ~ 0.8(1 - q)X 2 (fim)Bj_(G) \/ 1 - [i 2 . (10) 



Note that the perpendicular magnetic field consists of two 
mutually perpendicular components B^ = B p + B v , where 
Bp is radial ( in the plane of a disk) component, and B v is 
the azimuthal one. 

The mentioned above parameter C, according to Silant'ev 
(2005), is described by formula: 

C = (I - q)rP (S> 2 ) f B 

= 0.64(1 - g )r 1 (T) A 4 ( M m)(i?' 2 )/ B /3. (11) 

(T) 

Here, is the mean Thomson optical length of turbulent 
eddies, the values 5' and B' denote fluctuating components of 
corresponding values (5 = So + 5' , B = B + B'). In Eqs.(9) 
and (10) the values <5||,(5_i_,B|| and B^ are the mean values 
of corresponding quantities, where, for brevity, we omitted 
the subscript "zero". The statistical averages (5') = and 
(B') = 0. The numerical coefficient fs — 1 is connected with 
the integral from corresponding correlation function of values 
B 1 in neighboring points of turbulent atmosphere. 

The observed degree of the light polarization and the posi- 
tion angle are derived from parameters (9) by the usual way. 
For particular cases of pure normal (5± = 0) and pure per- 
pendicular (<5|| = 0) magnetic fields expressions (9) can be 
derived analytically. For the first case we have: 
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It is interesting that position angle x weakly depends on pa- 
rameter k/i < 1, which characterises true absorbtion in the 
atmosphere. For C>1 this term can be omitted. In the first 
equation this term can be also omitted, but the nominator 
p(//, 0) = p(/z, q, 0) strictly depends on degree of absorbtion q 
(see Silant'ev 2002). 
For perpendicular magnetic field are there the formulae: 
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For pure perpendicular magnetic field the position angle \ = 
is due to the axial symmetry of the problem. So, the electric 
wave oscillations in this case occur parallel to the surface of 
an accretion disk. The case x ^ can be only realized if there 
exists some Bu component. 

The degree of linear polarization p and position angle \ 
depend on parameters a and b in a fairly complex form. Of 
course, the magnitude of polarization decreases with the in- 
crease of a and b. The numerical calculations demonstrate 
that the relative polarization degree p(n, B)/p(/i, 0) is sym- 
metric function of parameters a and b. The position angles 
X do not possess this symmetry. The detailed description 
of parameters p and %, and their wavelength dependence is 
presented in Silant'ev et al. (2009). Note that the radial B p 
and azimuthal B^ components of the perpendicular magnetic 
field B^ integrally give the same degree of polarization. 

As it was mentioned, in Eq. (1) we have Bbh = #||- In 
this case we can to use more simple expressions following from 
Eqs.(12) : 
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If we know the inclination angle i, the observed polarization 
degree p(n, B) and the position angle x, then the expressions 
(14) allow us to know <5||, i.e. the value B\\, and the value 
1 + C — kfi for models of accretion disks with different pos- 
sible values of degree of absorbtion q. Note, that increase of 
Pr{^,q) occurs more fast than decrease (1 — q) with grow of 
parameter q. 

3 The correlation between Bbh and 
Mbh-, and optical polarization of 
AGNs and QSOs 

Using formulae (9), we calculated the polarization degree 
p(n, B) and position angle x( n i B) for the typical parameters 
of AGNs considered by Zhang et al. (2005). Our calcula- 
tions assume Mbh = 1O 8 M , inclination angle i = 60°, and 
effective band width with A res t = 0.55/xm. 

The first important step is to derive the scale length 
R\ defined by the radius in the accretion disk where the 
disk temperature matches the wavelength of the monitor- 
ing band. The standard accretion disk (Shakura and Sun- 
yaev 1973) is characterized by effective temperature profile 
T e = Tbh(Rh /r) 3 / 4 , where Rr is the black hole horizon ra- 
dius. There are the series of papers (Koshanek et al. 2006; 
Poindexter et al. 2008; Morgan et al. 2008), where the semi- 
empirical method of determination of the accretion disk scale 
has been developed. The authors used microlensing variabil- 
ity observed for gravitationally lensed quasars to find the ac- 
cretion disk size, and observed (in rest frame) wavelength 
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relation. It is very important that the scaling is appeared to 
be consistent with expectation from the geometrically thin 
accretion disk model of Shakura and Sunyaev (1973). 

Poindexter et al. (2008) have presented the following scal- 
ing size- wavelength relation: 

/ \ \4/3 

Rxicm) = 0.97- 10 10 x 

\lfxm ) 

(M BH \ 2/3 ( L bol \ 1/3 

Here R\ is the radius in accretion disk that corresponds to 
the observed effective wavelength A, L bo i is the bolometric 
luminosity and L E dd = 1-3 • 10 38 (M BH /M Q ) erg s 1 is the 
Eddington luminosity, e is the rest-mass radiation conversion 
efficiency The commonly accepted relation between the Lbol 
of the accretion disk and the accretion rate M is Lbol = £ Mc 2 . 

First of all we estimate the depolarization parameter a from 
Eq. (10) (the parameter b, in this approximation B(R\) = 
B\\(R\), is equal to zero). We calculate a suggesting the 
power-law radial dependence of the magnetic field 

B(R x )=B BH (^y . (16) 

The radius of the black hole horizon is (Novikov and Thorne 
1973): 

R H (cm) = ^£(1 + V / l 3 ^!), (17) 

where a* is the spin of the Kerr black hole. 

For our calculations we use the fundamental relation be- 
tween the bolometric luminosity and the mass of a super- 
massive black hole in AGN (see, for example, McGill et al. 
2008; Ziolkowski 2008, Denney et al. 2009). According to 
Ziolkowski (2008) we have 

M BH (g) = 5.71 • 10 7 L° o f 5±O ' O36 (44)M , (18) 

where Lb ;(44) = L^/IO 44 ergs -1 . Let us remind also the 
very important relation Lbol = 9-AL(5100X), where L(5100X) 
is differential luminosity at A = 5 100 A 

Eqs. (1), and (15)-(18) allow us to obtain the magnetic field 
B(R\) for power-law index n: 

B(R X ) = B {] (R X ) = 

10 9.26-2.036„-0.81 S+ 0.055 S n A -4n/3 /n( ^ )5 (jg) 

where s = \og{M B H /Mq) and 

/nM-^a + v/r^fr- (20) 

The analogous formulae can be derived for other models (see 
Eqs. (2)-(5)). Note that in all the formulae magnetic fields 
are taken in Gausses, the radii R B and R\ in cm, wavelengths 
in /zm, and mass in grams. To obtain the parameter 8\\ we 
substitute the formula (16) to first expression in Eq.(10). 

The coefficient /„(a*) is tabulated in Table 1 for a number 
of values of parameters e and a* . In this Table we used the 
known relation between the spin of the Kerr black hole a* 
and the accretion radiation efficiency e (see Tables in Krolik 
2007, and Shapiro 2007). 



Table 1: The values of the coefficient f n (a*). 
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Table 2: The values of the parameter 5^. 


a*; e 


0; 0.057 


0.5; 0.1 


1; 0.42 - 0.9; 0.032 
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= 0.5 


30.337 


32.182 


30.586 23.347 
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= 0.81 


11.463 


12.613 


11.615 7.500 
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= 1 


6.313 


7.104 


6.417 3.739 
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= 5/4 


2.880 


3.338 


2.939 1.496 
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= 3/2 


1.314 


1.568 


1.346 0.599 
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= 2 


0.273 


0.346 


0.282 0.096 



For the same values of parameters n,e and a*, as in Ta- 
ble 1, and the values X re st — 0.55/im, M BB — 1O 8 M , we 
calculated first the parameter <5|| (see Eqs. (10)). The results 
are presented in Table 2. The parameters 8\\ for other wave- 
lengths can be calculated from (5||(A = 0.55(Um) by the simple 
relation: 

/ \ \2-4n/3 

*||(A2) = *||(Ai)(^J • (21) 

The Faraday depolarization is not effective if the param- 
eters <5|| is small, 5ii < 1. From Eq.(l) we see that increase 
of ,s = M B h/Mq decreases magnetic field Bn, i.e. the de- 
polarization parameter 6\\ are to be also decrease. This is 
confirmed by numerical calculations. So, for s = 5 parameter 
5\\ lies in the interval (7000 - 40), depending on values n, a» 
and e. For s — 10 this parameter acquires the values between 
1 and 0.01. The largest value of parameter 5u corresponds 
to smallest value of n. The values of 5u , corresponding to 
n = 0.5, are to 100 - 200 times larger than those, correspond- 
ing to n = 2. It is clearly why this occurs. Large n give 
rise to small values of magnetic fields B(R\), as compared 
with those corresponding to small values of n. Dependence 
of <5|| on parameters a» and e is not so drastic as the s and n 
-dependencies (see Tabic 2). 

After that we assumed the inclination angle of an accretion 
disk i = 60° and have calculated the parameter a according to 
relations (10) for conservative and non-turbulent atmosphere 
(q = 0, C = 0). Using Eq. (12), we calculated the polarization 
degree p and position angle x- The results of calculations 
are presented in Table 3. Remember that position angle x 
is equal to zero when electric wave oscillations are parallel 
to the accretion disk plane. Remember also that for i — 
60° the polarization degree pt(h = 0.5) is equal to 2.252% 
(see Chandrasekhar 1950). Our calculations show that for 
relation (1) the effect of Faraday depolarization acts efficiently 
in the case of standard accretion disk for n < 2. For any 
distribution of magnetic field with n > 2 the depolarization 
is very low, and the polarization parameters correspond to 
the classic limit Pt(m) (Chandrasekhar 1950). This means 
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Table 3: The values of polarization degree p%, and position angle 
X° for inclination angle of an accretion disk i = 60° 
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0.68%; 36° 
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1.06; 31 
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1.16; 30 


1.27; 28 


1.80; 18 


n = 
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1.88%; 17° 


1.77; 19 


1.87; 17 


2.16; 8 


n = 


2 


2.23%; 4° 


2.22; 5 


2.23; 4 


2.25; 1 



that there is a direct relation between the boundary value of 
n (where depolarization does not yet occur) and the black 
hole mass magnitude. So, for Mbh < 10 8 M Q the strong 
depolarization exists for n < 1.5 — 2 (see Table 3). 

The technique of estimation of observed polarization from 
optically thick magnetized accretion disk, suggested in this 
section, uses relation (12). If the observed polarization dif- 
fers strongly from these estimations, it means that formula 
(12) does not valid, i.e. the assumptions used by Zhang et al. 
(2005) are also invalid. The same technique can be also used 
for other models of accretion disk structure. So, the main 
problem to use this technique is connected with the observa- 
tion of polarization degree p and position angle \ from AGNs. 
This is separate, fairly difficult, problem, especially for esti- 
mation of position angle and its orientation with respect 
to accretion disk plane. We do not discuss these problems. 

Strictly speaking, our calculations of radiation polariza- 
tion from the ring of the radius R\ can be used directly if 
the gravitational potential GMbh / R\ in the observed ring is 
much smaller then that near the vicinity of black hole. The 
strong gravitational fields near the black holes influences the 
Stokes parameters of outgoing radiation when it propagate 
to a distant observer. The detailed calculations of this effect 
have been made by Connors et al. (1980), and Dovciak et 
al. (2004). According to these papers, the results of calcu- 
lations of the Stokes parameters J, Q, and U the local (rest) 
reference frame in accretion disk are to be integrated along 
the geodesic paths of photons escaping from every point of 
the ring R\ . This integration gives rise to the deformation of 
the spectrum of radiation (reddening), and to the change in 
Stokes parameters due to the rotation of polarization plane. 
What is important is that the degree of polarization does not 
change its value. As a role, the distance R\ for optical wave- 
lengths is far from the black hole, and our consideration is 
valid without the gravitational corrections. 

Frequently the constitution of black hole vicinity is more 
complex than the existence of accretion disk, one can exists 
some types of radiating jets, toroidal clumpy rings etc. The 
considered mechanism of light polarization is only one from 
possible mechanisms. Besides, usually the observed position 
angles \ are not connected directly with the plane of accretion 
disk. In this case we are to check that the differences of 
observed position angles Xofcs(Ai) — Xobs(Xj) coincide with the 
corresponding theoretical differences for different wavelengths 
Xi and Aj.This procedure makes the estimation of magnetic 
field and unknown disk's inclination angle i more distinct then 
for the single wavelength observation. 



3.1 Application of theory to the galaxy 
NGC 4258 

Now, as an example of using our formulae, consider the Seifert 
II galaxy NGC 4258. It is found that this galaxy has ac- 
cretion disk (see , for example, Modjaz et al. 2005). The 
accretion disk has an almost edge on orientation with incli- 
nation angle i ~ 83°. In the center there is a black hole with 
M BH ^ 3.9 • 1O 7 M (see Herrnstein et al. 1999). The po- 
larized continuum radiation from the nucleus of NGC 4258 
was observed by Wilkes et al. (1995). We consider here only 
observation of continuum polarization at A ~ 0.55/im. The 
polarization degree is equal to p = 0.35 ±0.01%, and position 
angle relatively the disk's surface x — 7±1°. For conservative 
atmosphere the polarization px{i = 83°) = 6.9%. 

Small polarization degree shows that <5|| are to be great. In- 
deed, from first Eq.(12), taken for conservative non-turbulent 
atmosphere, we obtain <5|| ~ 161.56, that corresponds to 
B\\{R\) — 667. 6G and 5\\/j, = 19.69, which gives, according 
to second Eq.(12), the value x — 43.55°. This angle is far 
from observed value 7°. So, it is impossible to explain the 
observational data , assuming that accretion disk is the con- 
servative non-turbulent atmosphere. 

Now, we suppose that accretion disk is turbulent (C ^ 0) 
and conservative medium (q = 0). In this case Eqs.(14) give 
rise to solution: 6\\ = 39.135, and C = 18.129. Using Eq.(10) 

and (11), we obtain B\\{R X ) = 161. 7G, and t[ T) (B a ) ~ 929 
(for estimation we take /# = 1). We do not know the charac- 

(T) 

teristic dimension of turbulent eddies. If we take r{ =0.1, 
then we obtain B' = ^(B' 2 } = 96G. For t[ T) = 0.5 the value 
of B' will be to 5 times lesser: B' ~ 19G. As far as it goes 
the g-dependence we know that the absorbing atmosphere has 
more great polarization degree pt(p-) that the usual conserva- 
tive one (see Silant'ev 1980). The term fc/i in Eqs.(l) is small 
value (remember that k < 1 and fi = 0.122). So, according to 
first Eq.(14), the <5|| will be greater than that in conservative 
case. Neglecting small term kfi in the second Eq.(14), we see 
that the level of magnetic field fluctuations C also increases 
its value. 

It is important that our value B^(R\) ~ 161. 7G we ob- 
tained directly from observational data and our asymptot- 
ical formulae (14). We do not use the basic relation (1). 
The value B\\(R\) from Eq. (19) was derived with the us- 
ing this relation. For our case s = 7.591, A = 0.55^m, and 
B»(R X ) = 162G, Eq. (19) acquires the form: 



fn — 



B(R x ,0) 



0.1254- 10 1 - 2723 ™. 



(22) 



For conservative atmosphere and n = 0.5 we have /0.5 ~ 0.54. 
For n = 0.81 one has /o.si — 1-34, for n = 1 the corresponding 
value /„ ~ 2.35. Further, with the increasing n the values 
/„ increase exponentially. This behavior is opposite to that 
presented in Table 1. According to Table 1 and definition 
(20), every values of /„ are equal to 1 at n = 0, and then they 
diminish monotonically with increasing of n. Our values of 
/„ in Eq.(22), obtained from observational polarization and 
relation (1), are equal to 0.1254 at n = , and then they 
increase very rapidly with increasing n. 

This means that for NGC 4258 relation (1) distinctly does 
not takes place. True absorption increases 5\\ and, corre- 
spondingly, the value B(R\,q) > B(R\,0). The absorbing 
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atmosphere increases the /„-values in Eq.(22). 

Note that, according to selfconsistent model of Pariev et 
al. (2003), parameters n < 1 correspond to gas density of 
accretion dusk which increase with the distance from black 
hole. It seems fairly unphysically. For n > 1, the discrepancy 
between theoretical /„ from Table 1, and the values (22) is 
most large. It means that it is necessary to take into account 
the turbulent magnetic field into accretion disk and its con- 
tribution can radically change the relation between magnetic 
field strength at the horizon of a massive black hole and its 
mass determined by Zhang et al. (2005). 

4 Wavelength dependence of polar- 
ization 

The most important information on the magnetic field distri- 
bution in an accretion disk can be obtained from the wave- 
length (frequency) dependence of polarization degree and po- 
sition angle. This dependence follows from Eqs.(12) for the 
case of pure vertical magnetic field B\\, and from Eqs.(13) 
- for the case of perpendicular magnetic field B±. The ex- 
pressions for polarization degree in both cases are analogous, 
but the formulae for position angles are different. According 
to second Eq.(13) for pure perpendicular magnetic field the 
spectrum of x(A) 1S almost flat. This is can be used, to some 
extent, for qualitative characterizing is there great perpendic- 
ular magnetic field, as compared with B\\, or not. In general 
cases we are to use the numerical calculation of (Q) and (U) 
from Eqs.(9). Because absorption degree q depends on wave- 
length, it is difficult to speak about wavelength dependence 
of p(X) and x(A) in general. For this reason, we restrict our- 
selves by the case of conservative atmospheres, where q = 
and k — 0. Such atmospheres are frequently used in analysis 
of polarization data. 

Remember that parameters S\\,± are proportional to A 2 , 
and the turbulent parameter C is proportional to A 4 . The 
observed wavelength, according to Wien's displacement law 
A = Const /T e , is connected with effective temperature T e . 
In most used accretion disk model of Shakura & Sunyaev 
(1973) T e - r~ 3 / 4 . In others models one generally considers 
T e ~ r~ s , where r is distance from black hole. According 
to these relations, we have R\ <~ A 1 ^, where R\ is distance 
corresponding to emission of A-radiation. Using Eq.(16) we 
conclude that B(R X ) ~ \~ n/s . If %_l > C > 1 the po- 
larization degree spectrum is characterized by the following 
expression 



p(n,B) 



1 



'11,-L 



\2-n/s 



(23) 



For most used value s = 3/4 and n = 5/4 (see Pariev et al. 
2003), we obtain p(A) - A" 1 / 3 . If one exists C > 5\\ >± > 1, 
the p(X) - spectrum is described by another formula: 



1 



p(n,B) - - 



1 



C X A {B' 2 )' 



(24) 



We do not know how (B' 2 ) depends on R\. One can consider 
two limiting cases. If (B' 2 ) = Const, then p(X) ~ A~ 4 . If 
(B' 2 ) ~ B 2 , then p(X) ~ A- 2 < 2 -"/ s ), i.e. this is square of 
spectrum (23). In any case the most sharp dependence, as 



compared with Eq.(23), demonstrates that there exists mag- 
netic turbulence. 

Now let us say some words about the spectrum of position 
angle X- For C » 1 the second equation (12) acquires the 
form 



tan2% 



C 



X -(2+n/s) 

{B' 2 ) ' 



(25) 



For (B /2 ) = Const we obtain from Eq.(25) x(A) ~ 
arctanA _ ( 2+ "/ ;s ). This gives for n = 5/4 and s = 3/4 the 
expression x(A) ~ arctan A~ n / 3 . For second limiting case 
B' ~ B\\ we have x(A) - arctan \-(^/^). 

Qualitatively the behavior of position angle x can be un- 
derstood from Eq.(12). If 8\\ » C > 1, we have tan 2% > 1, 
and x — * 45°. In this case x(A) does not practically depend 
on A. The existence of turbulent extinction C diminishes the 
position angle. In the limiting case x 1, we can use the 
known relation tan2x — 2x- In this case Eq.(25) directly 
presents the spectrum of x(A). So, for (B 12 ) ~ Const, the 
position angle x — > very rapidly, as ~ A~ n / 3 for n = 5/4 
and s = 3/4. For B' ~ B\\ the position angle tends to zero 
slowly, as ~ A -1 / 3 for the same n and s. Note that spectra 
p(X) and x(A) help us to estimate the inclination angle i, as 
this was explained in paragraph 3. 

5 Conclusions 

We have showed that the magnetic field strength - the black 
hole mass correlation can be probed by the optical polari- 
metric observations. The basic idea of this probe is taking 
into account Faraday rotation of polarization plane of radi- 
ation scattered by free electrons in optically thick accretion 
disk. Faraday rotation gives rise to partial depolarization of 
outgoing integral radiation, and to arising of characteristic 
wavelength spectra of polarization degree and position angle. 
Both effects are described by two depolarization parameters 
a and b (see Eqs. (10)) related with normal B\\ and tangen- 
tial B± components of magnetic field inside an accretion disk. 
The Faraday effect is stronger if the parameters a and b are 
greater than unity. In this case there are the simple relation 
between the power-law index n of magnetic field distribution 
inside the accretion disk and black hole mass. As a rule, 
the greater the central black hole mass Mbh the lower the 
power-law index. 

For massive Mbh ~ (10 8 — 10 9 )M© one should expect the 
polarization degree corresponding to classic electron scatter- 
ing without noticeable Faraday depolarization. Only for ex- 
tremely smooth (n < 1) magnetic field distribution one can 
expect a display of depolarization effects. For low massive 
M B h ~ (10 5 - 10 6 ) M black holes the depolarization effects 
can be displayed even for steep distribution of magnetic field 
n > 1.5- 2. 

As an example of application of presented theory we consid- 
ered the AGN NGC 4258, where is found the accretion disk, 
which certainly plays the main role in polarization emission 
in continuum. It was found that for this source the basic re- 
lation (1) between magnetic field and the mass of black hole 
does not takes place. This negative result gives rise to basic 
question for what sources the correlation relation (1) really 
exists and why. 
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